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Summary 


The  equations  of  motion,  for  lateral-directional  dynamic  stability  and 
response  of  tandem  helicopters  in  forward  flight,  are  formulated.  Theoretical 
expressions  are  developed  for  the  various  stability  derivatives,  and  the 
results  of  preliminary  calculations  for  an  HUP  type  aircraft  are  discussed. 


Introduction 


This  report  is  the  result  of  preliminary  studies  of  lateral -directional 
dynamic  stability  characteristics  of  tandem  helicopters  in  forward,  flight  *  The 
work  is  part  of  the  third  phase  of  a  program  of  theory  and  static  and  dynamic 
flight  testing  to  develop  a  practical  working  method  of  stability  estimation 
for  tandem  helicopters.  The  first  two  phases,  involving  longitudinal  stability 
are  essentially  complete  (Ref.  1  and  2).  The  work  reported  here  is  to  be 
followed  by  analysis  of  static  lateral  flight  tests  to  verify  the  theoretical 
lateral,  derivatives*  and  then  by  lateral  dynamic  response  tests  for  correlation 
of  those  characteristics  with  theoretical  predictions. 


Axis  Systems  and  Notation 


The  equations  of  dynamic  motion  developed  herein  are  in  terms  of  forces 
and  moments  along  and  about  body  axes  fixed  in  the  aircraft,  and  rotating  and 
translating  with  it.  Their  origin  is  at  the  aircraft  C.&.  •  the  X  axis  lies 
in  the  plane  of  symmetry,  perpendicular  to  the  rotor  shafts;  and  the  2T  axis 
is  parallel  to  the  drive  shafts. 

It  is  convenient  to  use  existing  expressions  for  the  rotor  forces  and 
moments  in  terms  of  axes  set  with  the  local  wind  seen  by  the  rotors .  Quantities 
measured  in  this  "wind"  axis  system  are  subscripted  (  )w  :  such  as  Tw, 
acting  along  the  axis  of  No-Feathering  (ANF);  and  the  side  drag  forces,  Yvv 
and  Hw  ,  acting  in  directions  perpendicular  and  parallel  to  the  plane  containing 
ANF  and  the  instantaneous  relative  wind. 

In  relating  the  two  above  systems,  it  is  convenient  to  use  an  intermediate 
system  subscripted  (  )i  ,  in  which  Tl  is  along  the  ANF,  and  Yi  and  Hl  are 
perpendicular  and  parallel  to  the  plane  of  the  ANF  and  the  X  body  axiB. 

In  the  body  axis  system,  of  course,  T  acts  along  the  rotor  drive  shafts, 
and  Y  and  H  are  perpendicular  and  parallel  to  the  X  axis . 

Symbols  used  in  the  theoretical  analysis  are  detailed  below: 


General 

(  )' 


Primed  quantities  refer  to  the  rear  rotor. 


±(  ) 
i(  ) 

A(  ) 


The  upper  of  the  two  signs  preceding  a  quantity  refers  to  the  front 
and  the  lower  to  the  rear  rotor  value. 

Algebraic  sura  of  the  front  and  rear  rotor  contributions  of  a  given 
quantity. 

Difference  between  fromt  and  rear  rotor  contributions  of  a  given 
quantity,  the  rear  rotor  contribution  being  subtracted  from  the 
front . 

Root  of  characteristic  equation. 

Coefficients  of  characteristic  equation. 


R  Routh's  Discriminant  for  characteristic  equation. 

C  )  Bar  denotes  ^  (  )  . 

(  )^  *  (  )_jus  Subnotation  denotes  fuselage. 

(  Subscript,  following  standard  practice,  denotes  partial  different ia- 

^  tion  with  respect  to  subscript  quantity;  (  . 

( — y - denotes  differentiation  with  respect  to  time;  i”’  )® 


D(  ) 


Derivative  with  respect  to  non-dimensional  time',  D(  ) 


:c  y: 


3- 


Forces  and.  Moments 


T 

H 

Y 

L,N 


Rotor  thrust,  or,  upward  component  of  rotor  resultant  force  (lbs.). 

Horizontal  component  of  rotor  resultant  force,  directed  aft  (lbs.). 

Side  component  of  rotor  resultant  force,  positive  to  the  right  (lbs.). 

Rolling  and  yawing  moments  of  helicopter  about  C.G.;  roll  is  positive 
for  right  side  down  and  yaw  for  nose  right  ( ft. -lbs.). 


Thrust  coefficient. 


Side  force  coefficient. 


Cn 

(Cy)t 


Horizontal  force  coefficient.  CH  -  J^TTR^C SK)* 

Side  force  coefficient.  Cy 

Rolling  moment  coefficient.  “  C-as?)* JL 

Yawing  moment  coefficient.  0.^  — 


Fuselage  side  force  coefficient.  CO  Y  j  r  yp\ 

Fuselage  rolling  moment  coefficient.  CQ-jl  R-1  V"1  jL 

Fuselage  yawing  moment  coefficient. 


Physical  Characteristics 

Id  Number  of  blades  per  rotor. 

C  Average  blade  chord  ( ft . ) . 

C  Flapping  hinge  offset,  distance  from  blade  flapping  hinge  axis  to 

rotor  shaft  centerline  (ft.). 

)c\0  Average  height  of  front  and  rear  rotors,  measured  from  fuselage 

longitudinal  axis  (ft.). 

TtSkuSla.  Mass  moment  of  inertia  of  blade  about  the  flapping  hinge  (slug-ft?) 

T  T  T  Ta^  =  J.  Ara 

la,!}  Mass  moments  of  inertia  of  helicopter  about  its  x,  y,  and  z  body  axes 
0  respectively. 

Mass  product  of  inertia  of  helicopter  for  the  x  and  z  body  axes. 

^  OAT 

^  Radius  of  gyration  of  helicopter  about  its  (  )  body  axis. 

Linkage  term  representing  radians  of  swashplate  tilt  per  inch  of 
lateral  stick  motion. 


Linkage  term  representing  radians  of  swashplate  tilt  per  inch  of 
rudder  pedal  displacement . 


4. 


SI 

L 

~TfV 

Ms 

r 

R 

St 

W 


Distance  between  rotor  shafts 
Average  distance  from  C.G.  to 
Mass  of  helicopter  (slugs). 

Static  mass  moment  of  blade  about  the  flapping  hinge  ( slug -ft.) 

Ms  *  £  f  dm 

Radius  of  blade  station  or  element. 

Blade  radius  ( f t . ) . 

Area  of  vertical  tail  surface . 

Weight  of  helicopter  (lbs.). 


( f  t . ) . 

the  rotor  shafts 


•  1-  ki 


(ft.) 


Velocities  and  Angles 


V 

II 


Relative  wind  velocity  (ft. /sec.). 

Velocity  of  rotation  of  rotor,  clockwise  rotation  for  front  rotor 
and  counterclockwise  for  rear.  (rad. /sec.) 


o<  Angle  of  attack;  defined  for  body  axis  system  as  the  angle  between 

the  relative  wind  and  a  normal  to  the  rotor  shaft,  and  for  wind 
axis  system,  as  between  the  ’  elative  wind  and  a  normal  to  the  ANF. 

(3(  )  Sideslip  angle  of  helicopter,  positive  for  slip  to  right  (  •»!.). 

Roll,  positive  for  right  side  down,  and  yaw,  positive  for  nose  to 
right,  angles  of  helicopter,  respectively  (rad.). 


0*1 


Non-dimensionalized  roll  and  yaw  rates,  respectively. 

fepr  ^ “  b?- 

Longitudinal  tilt  of  tip  path  plane,  positive  aft  for  both  front 
and  rear  rotors  (rad.). 

Lateral  tilt  of  tip  path  plane,  positive  to  left  for  front  rotor 
and  to  right  for  rear  (rad.). 

Coning  angle  (rad.). 

Lateral  cyclic  pitch  of  front  rotor,  positive  to  right  (rad.). 
Lateral  cyclic  pitch  of  rear  rotor,  positive  to  right  (rad.). 


5S 

Bis 


Longitudinal  cyclic  pitch  of  front  rotor,  positive  forward  (rad.). 
Longitudinal  cyclic  pitch  of  rear  rotor,  positive  forward  (rad.). 


— Gfc —  —Co-llective-pitch  -anale^of-.-rotor-blade  s_(rad. ) .  _ _ 

(See  section  on  "Control  Effectiveness",  equations  13,  etc.) 
Incremental  lateral  cyclic  pitch  caused  by  positive  control  stick 
displacement  So. (rad. ) . 


I 


V 

I. 


I 

I 


Sx 

AiS  +  Kr  Sr 


Sr 


Lateral  control  stick  displacement,  positive  to  right  (in.)* 

(See  section  on  ’’Control  Effectiveness",  equations  l4,  etc.) 
Incremental  lateral  cyclic  pitch  caused  by  positive  rudder 
displacement  (rad.). 


Rudder  control  displacement,  positive  left  pedal  forward  (in.). 
Aerodynamic  and  Non-dimensional  Parameters 


/* 

A 

r 

a 

Q-t 

§ 

Ho 

K 

<r 

Y 

X, 

p * 

k 


il  n  V_ 

m 


Advance  ratio. 

Rotor  inflow  ratio. 

Air  density  (slugs /ft .3) 

Slope  of  blade  two-dimensional  lift  curve. 

Lift  curve  slope  of  vertical  tail  surface. 

Profile  drag  coefficient  of  rotor  blades. 

Induced  velocity  factor  to  account  for  the  effect  of  front  and 
rear  rotor  flow  interference. 

Aerodynamic  factor  to  account  for  the  change  in  lateral  flapping 
due  to  a  triangular  fore  and  aft  induced  velocity  distribution. 

*.  be 
TT  R 


Rotor  solidity.  (S' 

Locke1 6  blade  inertia  coefficient 


•  V 

4  I  BWE 

Blade  mass  constant.  v  -  Aac 


Non-dimensional  time  parameter 
Helicopter  density  parameter. 


Ms 

'V* 


P 


ft*  ~ 

Non-dimensional  roll  inertia  paraneter.  q 


as  yu* 

Non-dimensional  yaw  inertia  parameter.  | 


(ir 

(fr 


Non-dimensional  product  of  inertia  parameter. 


X 


_L  f  I 

ous  ' 


iB. 


6. 


Analysis 


Rotor  Forces  and.  Axis  Transformations 

The  rotor  forces  are  calculated  on  the  basis  of  several  simplifying 
assumptions,  as  follows: 

1.  The  induced  velocity  over  each  rotor  disc  is  constant.  The  inter¬ 
ference  induced  velocities  are  accounted  for  in  an  approximate 
manner,  based  on  longitudinal  stability  studies.  The  importance 
of  these  assumptions  is  less  than  in  the  longitudinal  case,  and 
they  are  believed  valid  where  only  perturbations  of  side  forces 
are  involved. 

2.  Assumptions  of  quasi-static  blade  motion  are  adopted.  The  blade 
dynamic  modes  of  motion  are  thereby  ignored,  except  for  the 
gyroscopic  lag  damping  effects.  The  effects  of  these  assumptions 
are  quite  generally  small,  since  the  modes  ignored  are  of  very 
high  frequency  compared  to  the  aircraft's  characteristic  modes. 

3.  Small-angle  assumptions  are  generally  adopted,  to  be  consistent 
with  the  small  perturbation,  linearized  equations  to  be  developed. 
Certain  higher  order terms  are  seen  in  the  development  to  be 
small,  and  are  ignored  in  the  interests  of  simplicity. 

The  basic  equations  for  the  rotor  forces,  in  the  wind  axis  system,  are 
given  below  (see  Ref.  3).  Where  double  signs  are  used,  the  upper  ones  apply 
for  the  front  rotor  (clockwise  rotation),  and  the  lower  for  the  rear  rotor 
(counter-clockwise  rotation). 

(cT)w  =  i(e„L  +  i(A)w 

-v  f  AU  +  (30  a,  -  Jju  A  l$o 


■+  ^  Jaaub,  -v-  y  ju*  Ab» 

w 


The  given  equations  (l)  reflect,  of  course,  the  induced  velocity  assumptions, 
and  the  ^£o  factor  in  (id)  allows  for  rotor  interference.  These  factors  must  he 
obtained  emperically  (HUP  interference  characteristics  are  contained  in  Ref.  l) 
and  they  represent  the  increase  in  rear  rotor  downwash  due  to  the  front  rotor  and 
the  smaller  upwash  component  on  the  front  due  to  the  rear  rotor* 


Similarly,  blade  flapping  is  given  by 
2a)  ( P>o)w  *  X  (  o  +  4  )\ 


I W 


2b) 


(a,) 


W 


|  §peo  -v 

1  JL  —  i-  hi 

4  %  r 


W 


2c) 

(bi)w  *  j 

3  (S  o  K 

i  +  -k  m* 

4- 


W 


V£© 

W 


$ 


Equations  (2)  of  course,  involve  the  quasi-static  flapping  assumptions,  and 
the  factor  K  in  equation  (2c)  allows  for  a  triangular  fore-aft  distribution  of 
induced  velocity  and  its  effect  on  lateral  flapping.  The  factor  K  is  given  by 

V  —  \  •+■  -2.  i5liw 

*  “  '  8  p  (So 

Equations  for  the  transformation  of  forces  between  the  axis  systems,  appro- 


priately  linearized. 

are 

3a) 

Ti  - 

Tw 

3b) 

X  - 

Yw  -  (3; 

3c) 

Hi- 

+  Yw  Pi 

and 

4a) 

T- 

fr\ 

II 

r1 

4b) 

Y- 

X  +  Tt  A,, 

4c) 

H- 

Hi  -  Ti  BiS 

Yw-  Hwft,  +  Tw  A, 5 


Mw  ■+■  Yw  S'L  "  Yw  6>3 


Aerodynamic  angle  transformations  are 


=  <=<*  -  ft 


l  *•  ft  4  # 


Blade  flapping  transformations  are,  similarly: 


L  Sic  —  (^l)w  BlS  "t 


I'/W  ~'t 


1)1  -  Ca,)w  -t  (Bi)w  ,3c 


(Bi)l  +  A  i 5  ^  +■  Als  -  (P-Ow  (3c 


(b»)t  =  (bi)w  ±  (a»)w  3c 


The  preceding  equations  are  sufficient  for  the  evaluation  of  the  derivatives 
of  the  rotor  forces  with  lateral  flight  variables.  These  derivatives  are  to  be 
-ta‘ken“atrHjhe~i“nl4ial™vaftues-“Oi™the— va-r-4'ab-le-S“>--™f©r-~e'X-ample»:™»----»---- — — - — — — . — — 


(  {U  ~  m<P  o  —  A  o  -  o 


9 


Before  proceeding  with  the  derivative 6 ,  however,  it  will  be  convenient 
to  express  the  inflow  ratio,  A  ,  in  terms  of  the  flight  variables  in  the 
body  axis  system.  Using  equations  (la)  and  (id) : 


M  tx  vv 


V  + 


m  c  j  > 
f  ( 


Substituting  (5a)  and  (5b)  into  (fa): 


Rotor  Side  Force  Derivatives 

For  the  side-force  derivative  with  sideslip,  from  (4b): 


aCy  a  (Cy)w 
d|3s 


(Cw)w  ■+•  A 


'S 


d  CCT )  < 

G>fc 


8a) 


turn 


Since  (Q)  w  is  a  function  of  Po  ,  ,  CX-i  ,  and  A  ,  which  in 

are  affected  by  /3s  ,  the  first  term  of  (6a)  can  be  expressed  as 


which,  from  (lb),  (2a),  (2b),  (2c),  and  (7b)  becomes 


~yM  A)S  F 

y  K  fco' 


where 


In  the  above  (8c)  and  (8d),  certain  higher  order  terms  have  been  dropped. 
Note  that  since  derivatives  are  to  be  evaluated  at  the  initial  condition  of 

symmetrtcttl--s-t-eady“-fdd•ght-~no”d^fstlTOtl~on7bf7^^^systgm71,^ 

necessary Tn' these  equations , 


The  second  term  of  (8a)  is  to  be  evaluated  'by  (lc) ,  which,  however,  may  be 
simplified  to 


-Lja  X  0O  -t-  J  \(Lt  •+  ^|xa,A 


The  last  term  of  (8a)  is 


Using  (8c)  and.  (8f)  in  (8a),  and  including  both  rotors, 


where  the  A  and  Z.  symbols  have  special  meaning  defined  in  the  Notation 
section,  and  f"  and  Ch  are  given  by  (8d)  and  (8e). 

It  may  be  noted  from  (8b)  that  the  derivative  of  side  force  with  rolling 
velocity  will  be  similar  to  the  derivative  with  sideslip,  except  by  a  factor 
and  the  addition  of  another  term  related  to  the  <J>  term  in  (2c);  It  may  there¬ 
fore  be  written  by  inspection  that 


It  is  apparent  in  (7b)  'that  the  derivative  of  side  force  with  yawing  velocity 
will  be  similar  to  the  derivative  with  sideslip,  except  by  a  factor  ,  and  a 


12. 


change  of  signs  for  the  rear  rotor.  By  inspection  of  (8g),  then 


-If  A  Ch 
V 


Rotor  Rolling  Moments 


Rolling  moments  due  to  the  rotor  are  of  two  kinds:  those  due  to  the  side 
forces  acting  at  J^0  above  the  C.G .  ,  and  the  hub  moments  associated  with  hinge 
Offset.  The  former  may  be  written  immediately  as 


10a)  _  ±  t  dO* 

df3s  A  X°  “5j£ 


10b) 


b<p 


_L  JL  c)C  y 

^  jto  c)  <fc 


5C  i  _  JL  JL  bCy 

"d  ^  JTo  d 


.  ...The  .hub. .moments  can  be’^rTtt'eh .  ly  aga^^dronlYrom^Referen-c'e'^rp-aR- 


11a) 


Q* 


> 


(90+  A)  +  (So  (  05  "  Iq  ^  <f>  ) 


... 

fa 


13 


The  derivatives  of  this,  using  equation  (6)  in  conjunction  with  (7) >  etc.  are 


11c) 


dCjL  = 

JL 

aCs> 

d'<f> 

V 

'bps 

1 


_  [£-) 

si  x&z  ''Jio  J 


Equations  of  Motion 

The  linearized  equations  of  motion  for  side  force,  rolling  and  yawing  moments, 
can  be  written  as  follows: 


+(  ICT  +  010  +  D)t  - 

15b>  (%  +^«D)Dt  +(^pD-^D^)<j>- 


(Cm. 
v  /U* 


A  0)D^V  +  (%  0  +  &5  DZ)4>  = 


18. 


and  o(Q  is  the  initial  angle  of  attack  of  the  X  axis. 

Expressions  for  the  aerodynamic  derivatives  in  equations  (15) ,  from  the 
previous  work,  are.  summarized  below; 


19 


22. 


18e)  -i\(r  i  cT 


Although  the  foregoing  equations  look  somewhat  formidable,  it  will  be  seen 
that  they  involve  many  factors  in  common,  so  that  their  computation  is  not,  in 
fact,  overly  lengthy.  It  is  to  be  noted  further,  in  the  next  section,  “that 
according  to  preliminary  calculations  at  least,  considerable  additional  simplifi¬ 
cation  is  permissible  for  typical  configurations. 


Sample  Calculations 

Preliminary  calculations,  for  the  HUP-1  helicopter  without  vertical  fins 
have  indicated  that  several  simplifications  of  the  above  equations  may,  in 
general,  be  permissible..  These  are  tentatively  listed  as  follows,: 


where 


l9J)  =  JLf*  +  +  f/'A3 

191^  o^,  >=  ^  p  G0  ■+  .sy*  A 
19n)  U  =  t  KyU  (3o 


A  few  preliminary  calculations  of  characteristic  modes  can  he  reported  at 
this  time.  The  helicopter  considered  was  the  HUP-1.  Tentative  results  are  as 
follows : 

l)  The  lateral  static  stability  is  positive  (  -f  dihedral  effect), 

<  ° 

g)  The  directional  static  stability.  ..is  negative  .  Oryg  <  °  .  due,  , to 

- — - uns-table-fiusela-g:e--ednt#ibUtion-.-~--! - - — 

3)  Damping  in  yaw  and  roll  are  positive  (  0v-r  and  0#p<O  ). 


4)  The  characteristic  modes  of  motion  for  these  derivatives  consist  of  two 
stable  real  roots  (simple  convergences),  and  an  unstable  oscillation. 
Both  real  roots  increase  with  increase  in  yM  .  The  period  of  the 
oscillation  increases  slightly  With  M  ,  whereas  the  damping  (cycles  to 
double  amplitude)  is  roughly  independent  of  . 

Studies  of  the  coefficients  of  the  characteristic  equation  indicate  that  for 
the  normal  HUP-1  configuration,  the  coefficients  can  be  simplified  to: 


Solutions  of  this  simplified  characteristic  equation  gave  approximately  the 
same  roots  as  the  more  complete  treatment  involving  much  more  complex  expressions . 


It  is  evident  _that  the  helicopter  can  not  be  dynamically  stable  with 
<  o  •  if  is  made  positive  (by  addition  of  vertical  tails, 


for 


example),  all  coefficients  above  in  (20)  will  become  positive,  and  the  stability 
will  be  determined 'by  Routh's  Discriminant, 
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K=  cs  c.c,  -  cucf  -  o.,cc* 


which  must  exceed  zero  for  stability.  With  approximate  coefficients  of  (20) 
however > 


21b)  R-  f  c. 


1-C,P 

1  J^5 


O 


\a  /  Cr 


V*/ 


in  which,  for^-in#  ^  O  ,  all  terms  contribute1  negatively.  This  denies 

hope  of  stabilizing  the  lateral  oscillation  merely  by  providing  static 

directional  stability. 

Solutions  of  the  more  complete  equation  have  shown  that,  as  On/3  is  made 
more  and  more  stable  (positive),  the  divergence  of  the  oscillation  is  reduced, 
but  not  overcome;  and  the  period  is  reduced. 

For  large  variations  in  the  other  derivatives,  the  approximate  equation  (20) 
is  not  valid.  More  complete  solutions  have  shown,  however,  the  following  effects: 

1)  Reduction  of  has  a  favorable  effect  on  the  unstable  oscillation, 

and  particularly  in  combination  with  positive  Oryj  ,  can  possibly 
stabilize  the  lateral  oscillation.  ' 

2)  Increase  of  damping  in  yaw  and  roll,  particularly  the  latter,  is 
beneficial . 

3)  Increase  of  side  force  derivative  is  beneficial,  though  possibly 

impractical .  ' 

Preliminary  conclusions  are  that  design  changes  to  improve  the  lateral 
oscillatory  instability  should  aim  to  a)  provide  static  directional  stability, 
b)  reduce  dihedral  effect,  c)  improve  damping  in  roll.  Further  calculations,  in 
greater  detail,  of  the  practical  possibilities  are  currently  being  undertaken,  in 
addition  to  experimental  verification  of  the  static  stability  derivatives  and 
dynamic  response  characteristics  discussed  above. 
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Conclusions 


Lateral  and  directional  static  and  dynamic  stability  derivatives  are  formu¬ 
lated  for  a  tandem  helicopter  in  forward  flight.  Equations  of  motion  are  written, 
from  which  the  dynamic,  stability  characteristics  are  predictable These  theoretical 
results  are  subject  to  experimental  verification >  and  should  be  regarded  as 
tentative . 


Sample  calculations  for  the  HUP-1  helicopter  indicate  negative  static  direct¬ 
ional  stability,  with  an  unstable  oscillatory  mode  of  lateral  motion.  It  appears 
that  provision  of  positive  static  stability  alone  can  not  stabilize  the  dynamic 
oscillatory  mode;  and  that,  in  addition,  reduced  dihedral  effect  and  increased 
damping  in  roll  are  required,. 
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